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Coherent dynamics of localized spins in semiconductors is limited by spectral diffusion arising 
from dipolar fluctuation of lattice nuclear spins. Here we extend the semiclassical theory of spectral 
diffusion for nuclear spins J = 1/2 to the high nuclear spins relevant to the III-V materials and 
show that applying successive qubit 7r-rotations at a rate approximately proportional to the nuclear 
spin quantum number squared (/ ) provides an efficient method for coherence enhancement. Hence 
robust coherent manipulation in the large spin environments characteristic of the HI-V compounds is 
possible without resorting to nuclear spin polarization, provided that the yr-pulses can be generated 
at intervals scaling as I~^. 
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I. INTRODUCTION 

The ability to enhance coherence times of localized 
electron or nuclear spin qubits in a semiconductor en- 
vironment is central for the development of a spin based 
solid state quantum computer^. Qubit coherence in these 
devices appears to be limited by the interaction with nu- 
clear spins composing the material-"'^, with other mech- 
anisms being less important. A detailed comparison of 
the many mechanisms contributing to localized spin de- 
coherence in semiconductors is given in Refs._35A Exper- 
imental evidence of the relative contributions of differ- 
ent decoherence channels is given for silicon in Refs.|ala. 
Here we consider qubit coherence under the strong effec- 
tive magnetic fields and low temperatures required for 
single spin read-out {B > 5 Tesla and T ~ 100 mK, see 
Ref. Q). Under these conditions, environmental nuclear 
spin fluctuations result primarily from inter-nuclear dipo- 
lar coupling, which becomes a source for time-dependent 
noise in the qubit Zeeman energy and for its consequent 
decoherence. This channel for phase relaxation is usually 
denoted spectral diffusion"'-^-*-^. 

It was pointed out recently that substantial nuclear 
polarization can suppress spectral diffusion^'^. However, 
the requirement for nuclear polarization adds significant 
device overhead, since polarizations on the order of 95% 
or more are needed in order to achieve sig nificant co- 
herence enhancement (See Fig. 3 of Ref. ll(Jl . Recently, 
several schemes for nuclear polarization were proposed^- , 
but it is still not clear whether these will be effective 
enough to compete with the nuclear spin diffusion rate 
[1 - 10 KHz, see Eq. ^ below]. 

Several pulse sequences have been designed to aver- 
age out the effect of nuclear-nuclear dipolar coupling, 
most notably WAHUHA and Lee-Goldburg schemesi^. 
Because these sequences are based on nuclear spin exci- 
tation, they inevitably generate a drastic time-dependent 
drift in the spin qubit Zeeman frequency. The amplitude 
of this drift can be as much as 50 G (typical hyperfine 
line- width for a localized electron in a III-V compound), 
fluctuating at the nuclear spin rotation frequency (typ- 



ically 10'* — 10^ Hz depending on the rf field intensity). 
The effect of these pulses is additional field fluctuation 
which adds to the decoherence rate and makes qubit 
control impossible. Under such a protocol, the fidelity 
of quantum operations would be drastically affected in 
a negative sense. Conversely, the nuclear spins caus- 
ing spectral diffusion are subject to an inhomogeneous 
field arising from the interaction with the qubit. This 
leads to a wide spectrum of nuclear resonance frequen- 
cies, requiring broad-band excitation in order to per- 
form nuclear spin rotation. Broad-band excitation usu- 
ally requires high power deposition, which is incompati- 
ble with the MilliKelvin temperatures needed for single 
spin measuremenlii^. These arguments suggest that any 
scheme for spectral diffusion control must rely on exci- 
tation of the qubits themselves, and not excitation of the 
nuclei. 

Here we consider the Carr-Purcell-Meiboom-Gill se- 
quence (CPMG)i^4i^ as a scheme for spectral diffusion 
control. The CPMG pulse sequence consists of the suc- 
cessive application of qubit 7r-rotations perpendicular to 
B at time instants {2n -\- 1)t, n = 0, 1, 2, . . .. This proce- 
dure leads to the formation of spin echoes at instants 
(2n -I- 2)t, whose coherence is substantially enhanced 
when compared to qubit evolution in the absence of tt- 
pulsing. Coherence enhancement takes place provided r 
remains below a threshold Tc, the typical time scale for 
environmental fluctuations. As r is decreased further, the 
coherence gain can be substantially amplified, with an in- 
creased cost in the number of n pulses. Here we develop 
a microscopic theory which is able to predict, without 
any fitting parameters, the inter-pulse time t„ needed to 
provide n times more coherence for a spin qubit subject 
to nuclear-induced spectral diffusion in a bath of general 
spin /. Our explicit calculations reveal that in realistic 
devices nuclear spin noise can only be suppressed if r is 
inversely proportional to the nuclear spin quantum num- 
ber / squared, establishing the overhead requirement for 
CPMG control of an isolated spin coupled to a nuclear 
bath of general spin /. 

Available spectral diffusion theories have focused solely 



on the case of / = 1/2 environments, using two-level 
telegraph noise models to analyze the coherence time of a 
central spin-l/3M. However many solid state spin-based 
proposals rely on high spin environments, an important 
example being the GaAs (/ = 3/2) quantum dot. Here 
one often has to add In (/ = 9/2) to enable single spin 
read-ouliiS^ in addition to a substantial amount of Al 
(/ = 5/2) to achieve confinement^^. InAs self assembled 
quantum dots are also promising with regard to optical 
manipulation and transpor^^iifiiiSiSii. The requirements 
for CPMG control in these structures are not yet known. 
In this work we first characterize spectral diffusion for 
/ > 1/2 environments and then analyse coherence control 
with CPMG. Our results estabhsh CPMG as a remedy 
for the absence of / = isotopes in III-V materials^ as 
long as existing schemes for fast and precise spin rotation 
are further develop ediS^i^iSi. 



gleet of flip-flop correlation a posteriori by noting that 
the condition Max Ft ^ 1 is never realized in the physi- 
cal cases considered here. 



A. Calculation of the flip-flop transition rates in a 
spin-/ bath 

In the semiclassical theory, nuclear spin fluctuation is 
decoupled from the qubit by setting 5*2 — > 1/2 in the full 
Hamiltonian Eq. (Q. This leads to an effective Hamilto- 
nian Ti for nuclear spin evolution under the inhonioge- 
neous field produced by the qubit. 



"H = "Hq 



^ nra \^)^ 



(2) 



II. SPECTRAL DIFFUSION DECAY INDUCED 
BY / > 1/2 NUCLEAR SPINS 



no = -Y.{iiB-\A 



Inz-^Yl ^" 



i^nz^rtiz ■) \'^} 



Our qubit/spin / bath Hamiltonian is given by 
n = jsBS,--fiBj2lnz 

n 

~r y^ -^n-^nz^ z ^ / ^ ^nm^nz^mz 
n n<m 

+ / ^ bnm\In+Im- + In-Im+)- (1) 



The first term describes the qubit two-level structure, 
7s and Sz being its gyromagnetic ratio and spin op- 
erator. The nuclear spin operators 7„ are coupled to 
the qubit through the parameters An. These are deter- 
mined by hyperfine coupling (for an electron spin qubit) 
or by dipolar interaction (for a nuclear spin qubit). The 
inter-nuclear dipolar coupling bnm is responsible for nu- 
clear spin fluctuations (for details see Ref . jj) . In partic- 
ular, the last term of Eq. IpJ induces "flip-flop" tran- 
sitions between pairs of nuclear spins. When such a 
transition takes place, the qubit level spacing shifts by 
2Anm = \An — Am\, leading to phase randomization and 
to decoherence. The dipolar coupled spin-/ lattice acts 
as a bath which is coupled to the qubit with strength 

In the semiclassical approach to spectral diffusion^-, the 
intricate nuclear spin evolution stemming from Eq. (^ 
is effectively described by a set of uncorrelated classi- 
cal stochastic transitions (flip-flops), with a characteris- 
tic rate F obtained from first principles. It is appropriate 
to use a sudden-jump model for the description of these 
processes because the time scale for thermal fluctuation 
is much shorter than i (fi/fcsT <C ns^ ^r ^ ms). Corre- 
lation between flip-flop events of fluctuating pairs located 
near each other become important only when these pairs 
flip-flop several times within r (or equivalently Fr 3> 1). 
For relevant time scales (2r < T2) we can justify our ne- 



In Eq. Q, we have introduced a flctitious frequency u 
which allows a connection to the method of moments^. 
After we calculate the flip-flop rate, we will take the limit 

Without loss of generality, we assume n — I, m — 
2. The eigenstates of Ho are given by |Mi,M2,...), 
with Mi = —I,—I+ 1,...,/. The flip-flop operators 
Fi2{t) do not connect all states because they conserve 
J' — (Ml -I- M2)/2. Hence flip-flop dynamics break the 
pair Hilbert space into disconnected subspaces labeled 
byj' = /,/— 1/2,...,— /. Each subspace has dimension 
2{I — \J'\) + 1 and can thus be alternatively labeled by a 
pseudo-spin /' = I — \J^\, which is seen to correspond 
to projection of the pseudo-spin operator (/i — /2)/2 
with projection index M' = (Mi - M2)/2 = -/', ...,/'. 
Each nuclear pair contains two subspaces for each /' = 
0, 1/2, 1, . . . , /— 1/2 (since J' can be positive or negative) 
and one subspace with /' = /. 

The Fermi Golden rule rate for a transition M' ^^ M' — 
1 between the state |M') = |Mi, M2, M3, . . . , M^} and 
\M' - 1) = |Mi - 1, M2 -h 1, M3, . . . , Mn) is given by 

F = 2^6?2 E PiM3,...,MN) 

M3,...,Mn 

X |(M' - l|/i_/2+|M')|' SiEl,,_, - El,, c.),(5) 

where p{M^, . . .) are Boltzmann probabilities for the un- 
perturbed energies El^,. It is convenient to define an 
auxiliary function 

piu)^ Y, p{M^,...,MN)5{Elt,_^-Elt,-Lo), 



A/3,...,J\/ 



3,...,ilJjV 



(6) 



which can be interpreted as a density of states. After 
evaluating the matrix element in Eq. Q we get 

r = 2TTbl2i2I - I' + M'){21 - I' - M' + 1) 

x{I' + M'){I'-M' + l)p{io), (7) 

where we recall that I' = I - \Mi + M2I/2 and M' = 
(Ml — M2)/2. We can now calculate explicitly the mo- 
ments of the function /9(a;). For example, J p{uj)duj — 1, 
while the first moment is given by 



^12 



ujp{uj)duj 



M3,...,Mn 



,Mn){E%_,-E^ 



M' 




r/A Regime (c) (Log) 



= -(A2-Ai)-4&i2(A/i-M2-l) 
+4 Y^ ibik-b2k){mk). 



(8) 



Here the thermal average (m/j) = dlnZ/dP' can be cal- 
culated from the partition function 



E 



I3'm 



(9) 



with /?' = n-yiB/kBT (for B :s> lT,n bnm/li - 10 G we 
can neglect the dipolar energy in the Boltzmann distri- 
butions). The second moment is given by 



'-12 



{lj — UJ12) p{uj)duo 
16 Y. (^^k - b2kf{{mk - (m,))2), (10) 



where {{mk — (mfe))^) = 9^ In Z/ 9^/3'. The fourth mo- 
ment g^2 can also be calculated, allowing us to show that 
the ratio 912/ (3^12) ~ 1- Thus the Gaussian function 
provides a reasonable fit to Eq. (0, 



p{uj) 



•\/27rK^, 



: exp 



{ll> - UJ12) 

2k^2 



(11) 



The final expression for the flip-flop rate for a nuclear 
spin pair (n, m) ia2^ 

Tm' = (21 - r + M'){21 - I' - M' + 1){I' + M') 



x(/'-M' + l) 



2TTbi 



■ exp 



ujT, 



9k-2 



(12) 



When /?' ^ 1, there is only a weak temperature and B 
field dependence in Eq. I|12|l . For example, T = 100 mK 
and _B = 10 T leads to /3' ^ 0.1. In this approximation 
we have {nik) ~ and {{irik — (mfe))^ ~ I{I + l)/3. This 
leads to the approximate expressions 



46„™(2M'-1), 



(13) 



FIG. 1: Behavior of 1/e decay times for the first CPMG echo 
in pseudo-spin subspaces I' < 5/2 as a function of flip-flop 
rate F. Circles: calculation with M'-dependent rates Tj^j,. 
Solid lines: calculation with average rate (F^')- Each curve 
achieves minimum coherence at F = Fc, which is identifled 
as a critical correlation rate. Small F (or strong coupling A) 
imply low frequency noise with time correlation [reginie(a)], 
while large F (or weak coupling A) imply motional narrow- 
ing or white noise behavior [regime (c)]. Note the strong 
I' dependence in the weak coupling region. Here we have 
Min (T2) Si 6/A(2/' + 1), and Fc ~ {21' + 1)^12. 



16 



1(1+ I) 



E 



(6. 



(14) 



Knm represents the amount of energy the dipolar system 
can supply for a nuclear flip-flop to take place. Nuclear 
spins located near the center of the electron wave func- 
tion are characterized by ZJnm ^ Unm , with exponentially 
small flip-flop rates. These are said to form a frozen core 
which does not contribute to spectral diffusion (a simi- 
lar effect has been reported in optical experiments, see 
Ref. l23f) . On the other hand nuclei satisfying tJnm ~ I'^nm 
possess appreciable flip-flop rates together with a sizable 
hyperfine shift A„m ■ These are located in a shell around 
the center of the electronic wave function, and are respon- 
sible for most of the spectral diffusion coherence decay^. 
Below we show that to a very good approximation the 
M' dependence in Eq. (|12|l averages out, so that we may 
consider just the average rate 



F(/', /) = (Fm'> « e(/', /) r(/, /) ^ /r 



/2 



(15) 



with the parameter ^(/',/) connecting different sub- 
spaces of a spin-/ pair. 



^{I',I) 



I' {I' + 1) [1 + 5/(2/ + 1) + /'(2/' - 10/ - 3)] 



/(/+l)[l + 2/(/+l)] 



(16) 



B. Random walk theory for the fluctuation of a /' 
subspace 

We now formulate a random walk theory for the cou- 
pling of the qubit to the flip-flop dynamics in each 
pseudo-spin subspace /' and then construct the coher- 
ence decay as a weighted average over all pseudo-spin 
subspaces. Our model assumes that each /' subspace 
fluctuates independently, i.e. pairs are uncorrelated. As 
a result, the qubit Zeeman frequency is subjected to sev- 
eral random fields of the form Act/' (t) , where ct/' is a 
hard wall random walk variable assuming the values 



CT/'(t) = {-2/',-2(/'-l),...,2/'}. 



(17) 



The theory for qubit dccoherence resulting from the two 
/' = 1/2 subspaces is described in Ref.|3 It has a simple 
structure because for /' = 1/2, ai'it) can be written as 
(— 1)^(*\ with N(t) a Poisson random variable with pa- 
rameter Tt. The generalization of this telegraph process 
to the multi- level noise present when /' > 1/2 is achieved 
by defining an occupation probability vector p[t) and its 
Markovian evolution operator U{t): 

pit) = iP-i',P-i'+i, ■■■,?!'), (18) 

= cxp {-tA) ■ p(0) (19) 

= U{t)-piO). (20) 

For first-order flip-flop transitions we have 

Pt = Ti pi^i - {Ti + Ti+i)pi + r,;+i Pi+i, (21) 



states CT. The n-th echo of the CPMG sequence can then 
be simply obtained from the n-th power of Eq. (|25|! "'"^'^"^. 
The difficulty of evaluating Eq. H25(l lies in ensuring 
a proper treatment of s(t). We propose here a matrix 
method that enables the incorporation of any step-wise 
constant s{t) explicitly. After expanding and rearranging 
the integrals we obtain 



VI' {2t) 



It 
^(^'\^ I dt'. 



dtk... / dti {[ait'j) ■ ■ ■ ait[)] 
Jo 

x[CT(tfe)---CT(ti)])^, 



dt[ 



(26) 



where r > ifc > ■ ■ ■ > ii > and 2t > t^ > . . . > f^ > r 
are partitions before and after the tt pulse. Using Eq. if^ 
and the Markovian identity U{t[ — th) = U {t'^ — t) -U [t — 
tk), Eq. f^ can be written as 



VI' {2t) 



1 



2/' + l 



(1, . . . , 1) • M_ • M. 



+ ■ 



(27) 



where the initial probabilities are simply p^r (0) — 1/{2I'+ 
1). The matrices M± are given by (after the substitution 






(28) 



Aij — —TiSi^ij + {Ti + r,;+i)(5ij- — Ti+iSi+ij. (22) 

We can determine the echo decay for all Markovian pro- 
cesses described by A using the time-ordered correlation 
function {tk > t^-i > • • • > ii) 

{ai'{tk) ■■■(TI' (ii)) = {I ■■■!)■ Hi' ■ U{tk - tk-i) 
■111' ■ U{tk-i - tk-2) • • • Sj- U{h) ■ p(0), (23) 

defined by correlations between values of the spin state 
matrix 

S,, - diag {-2/', -2(7' - 1), ... , 21'} (24) 

at times ti,t2, ■ ■ ■ ,tk- The Hahn echo is then calculated 
using a double average 



u//(2t) = / / exp i / s{t)(7r{t)dt 

\\ \ Jo / / <t,,(0)=<t/ ^ 

(25) 
where the echo function is s{t) = 1 for t < t and s{t) = 
— 1 for t > T (we set A = 1 for simplicity, since it will be 
recovered later as the unit of time). The inner average 
is over ct/' (t) trajectories starting from ct/' (0) = a. The 
outer average is a thermal average over all possible initial 



M+(r)=^(z)M dtkU{T-tk)^mk{tk), (29) 
where mi{ti) is obtained from the recurrence relation 



mi{ti)^ dti^i-Ei'^U{ti-ti^i)^mi^i{ti^i), (30) 
Jo 



mi(ti) = S7'-C/(ti) 



(31) 



Going to Laplace space the convolution integrals are sim- 
ply converted to multiplications. 



mi{r) 



e-''^mi{t)dt = Ej' ■U{r) 



(32) 



Using Eqs. (|28|l . (|29|l we obtain closed expressions for the 
Laplace transforms of the matrices M± , 



M_{r) = ^—^ . 

ri + iSi' ■U{r) 



M+{r) = U{r) 



1 - iSr ■ U{r) 



(33) 



(34) 



These general expressions for multi-level Markovian echo 
can be further simplified when written as a function of 
the eigenvectors and eigenvalues of A. 

Solving Eqs. ||2ZJ), §^-§^, using ^-^, allows 
evaluation of the dccohercnce time T2 from the 1/e decay 
of the Hahn echo, Eq. (|27|1 . An analytical expression for 
the Hahn echo is obtained using Mathematicc?^ . Fig. Q] 
shows T2 as a function of average flip-flop rate T for sev- 
eral /' subspaces. First we note that for all /', averaging 
over the M' dependences in Tm' is an excellent approxi- 
mation [resulting from the outer average in Eq. 1)25(1 ]. For 
each /', we identify a critical bath correlation rate Fc for 
which the minimum decay time T2 is found. This leads 
to a natural division into three different regimes which 
can be seen to possess different capabilities for coherence 
enhancement by CPMG sequences. We recall that ap- 
plication of a CPMG sequence to Hahn echo given by 
exp [— (2r/T2)''] transforms this to exp (— 2nr/TcpMG) 
where Tcpmg = \T2/{'2'T)Y~^T2 provides a measurement 
of the coherence enhancement which is obtained for all 
d > 1 (no enhancement results for d = l)i^. 

We refer to Fig. ^ for an illustration of regimes (a), 
(b), (c). Regime (a) holds for F/A < Fc. This is the 
regime of sudden-jump spectral diffusion (strong coupling 
A) which has oscillations in the echo envelope due to 
precession under a finite A. Assuming that at most two 
flip-flops take place, we obtain an expansion up to second 
order in (F/A), 



vp{2t) 



27^r+^^'+ a'""' 



X [1 + (2/' - 1)0-"] + ( ^ ) [l-cos^ 



21' - 1 



-e-''(l-cos25) 



(35) 



where p = 2tF and 5 = 2rA. When r < 1/A, CPMG 
provides coherence enhancement proportional to (r2/r)^. 
Regime (b) holds for F/A ^ Fc. This is the regime of 
critical spectral diffusion, which occurs when T2 reaches a 
minimum. The minimum is a result of the inevitable mix- 
ing of the random walk noise in a finite walk, such as that 
in the pseudo-spin subspace considered here. Coherence 
enhancement is approximately proportional to (T2/T) for 
r < 1/A. Finally, regime (c) satisfies F/A > Fc. This 
corresponds to the continuum or weak coupling regime. 
For 2r ^ T2 nuclear spin fluctuation is so fast that the 
coherence of the central spin-1/2 is enhanced. Motional 
narrowing takes place when the stochastic process tra- 
verses the random walk space several times; hence it is 
insensitive to all eigenvalues of A except for the low- 
est non-zero eigenvalue a* and associated eigenvector v* . 
The latter observation allows us to derive an analytical 
expression for the decoherence time at motional narrow- 
ing, when qubit spin echo vp w exp (— 2r/r2) with 
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2/' + l 
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FIG. 2: Coherence times T2 and inter-pulse times rioo re- 
quired for 100 X coherence enhancement as a function of nu- 
clear spin quantum number I. We find that both T2 and rioo 
scale as ~ 1//^'^. 



When Eq. l|l^ holds, no CPMG enhancement is ob- 
served. We note however that as r decreases below a 
certain threshold value r^ « l/(a*F) « {21' + 1)2/(8F), 
crossover to vi' ~ exp(— t"^) behavior takes place and 
consequently the CPMG pulse sequence will again pro- 
vide effective enhancement of coherence. The threshold 
value Tc is related to continuous Gaussian noise, since 
the limit F ^ 00, A — > and I' —> 00 with Tc and /'A 
constant is given by Eq. (All) of Ref. |3- 

In order to determine the dependence of spin coher- 
ence on the bath nuclear spin value J, we make an aver- 
age of Eq. (|27|l over all /' subspaces and nuclear species 
s. We recall that most materials considered here are a 
mixture of different isotopes. In a strong external fleld, 
only flip-flop between the same nuclear species needs to 
be considered, since flip-flop between different species is 
strongly suppressed due to an offset in hyperfine shifts. 
The resulting echo envelope is 



(t,),(2r) = ^p(5,/')t'/'[2r,C(/',/)r], 



(37) 



21' 



«i*E«:(S/')«A2 



oc/' 



(36) 



where p(s,/') oc f^ {21' + l)cosh[2{I - r)h'yiB/kBT] 
with fs the isotopic abundance for nuclear species s. 
Eq. (|37l) shows that substantial nuclear spin polariza- 
tion (or extremely low nuclear spin temperature T) would 
be required to suppress spectral diffusion. For example, 
p(0) = 0.9 (corresponding to 90% of nuclear spins hav- 
ing maximum magnetization /) increases T2 by a fac- 
tor of 2 only, while for the completely unpolarized case, 
p{I') = 2(2/'-|-l)/(2/-hl)2^ Numerical evaluation shows 
that the dependence of Eq. (|37ll on F/A is qualitatively 
similar to that seen for the individual /' subspaces in 
Fig. n but with different scaling behavior for the char- 
acteristic quantities. For example, Fc « [{21 + l)/3]^'^, 
allowing possible crossover between the three regimes as 



/ increases. 



III. NUMERICAL RESULTS 

The complete coherence decay is obtained as a prod- 
uct of decays from muhiple fluctuator pairs n, m. Most 
/ > 1/2 materials of interest have similar dipolar cou- 
plings {bnm) SO that we may conclude that CPMG con- 
trol depends primarily on the shape of the qubit wave 
function (which determines the distribution of A„m) and 
on the magnitude of the nuclear spin /. Fig. El shows 
our numerical results for T2 and tioo (inter-pulse time 
needed for 100 x coherence enhancement) for the case of 
donor impurities (electron and nuclear spin qubit) and 
a III-V quantum dot (electron spin qubit - see Ref. Q)- 
For an electron spin it turns out most nuclear pairs are 
fluctuating in regime (a) with TT2 <^ 1, implying the 
threshold for coherence enhancement Tc is quite close to 
T2. As T decreases below Tc ~ T2, coherence enhance- 
ment scales as (T2/r)^, yielding tiqo ^ 0.05r2. On the 
other hand, coherence decay of nuclear spin qubits is sen- 
sitive to fiuctuators in regimes (a)-(c) with rT2 < 1, 
which makes the threshold Tc somewhat smaller than 
Tj^. Here rioo ~ O.O3T2, showing a smooth crossover 
from 1/J to 1/P. Our study shows that the application 
of TT- rotations every r ~ 0.1 /its for an electron spin qubit 
and T ~ 10 /is for a nuclear spin qubit is sufficient to in- 
crease coherence by two orders of magnitude, effectively 
reaching the phonon emission limi^. 



IV. DISCUSSION AND CONCLUSION 

Fig. answers the previously open question of the de- 
pendence of spectral diffusion with increasing /. There 
had been earlier speculationSl that the spectral diffu- 
sion rate I/T2 might decrease with increasing I because 
of motional narrowing [due to faster flip-flops - notice 
that averaging over the /' subspaces in Eq. ((T3|) results 
in r (X I^] . We show here (see Fig. ^ that the drastic 
increase in size of the Hilbert space for increasing / does 
actually compensate for the strong amplification of flip- 
flop rates and that therefore the system does not enter 
the motional narrowing regime, resulting instead in more 
spectral diffusion induced decoherence as / increases. 



Another interesting point to note is that the repetition 
rate needed to suppress spectral diffusion (of the order of 
J^ X MHz) is significantly different than that correspond- 
ing to the strength of dipole-dipole interactions, which 
are given by the root mean square of the last two terms 
of Eq. n This root mean square average is commonly 
referred to as the square root of the second moment (see, 
e.g., Eq. 3.53 of Ref. Hal . For the systems of interest 
this is of the order of /xKHz. The rate for control- 
ling spectral diffusion is significantly faster than dipole- 
dipole coupling because hyperfine coupling (absent from 
the second moment equation) is very strong (up to MHz) 
and the spectral diffusion decay occurs due to the collec- 
tive fluctuation of several nuclei under the field of the 
electron. We emphasize that the relative time scale for 
CPMG control of spectral diffusion and its dependence 
on the nuclear spin value / that we have established here 
can only be derived by carefully considering the combined 
effect of hyperfine interaction and collective dipolar fiuc- 
tuation of IC* nuclear spins under the inhomogeneous 
hyperfine field produced by the electron, as we have done 
in this paper. 

The effectiveness of the CPMG sequence in controlling 
spectral diffusion depends largely on the availability of 
fast methods for precise spin manipulation. The pulsing 
time has to be much smaller than r, while the spin rota- 
tion angle should be precisely tuned to tt^*. The train of 
TT pulses must also be applied without heating the sample 
above the required ~ 100 mK temperature. Currently, 
there are several proposals for overcoming this techni- 
cal challenge, including all electrical spin resonancei^, 
optical manipulation^^, and implanting a low power mi- 
crowave source nearby the qubit^. If this technical prob- 
lem can be solved, our work indicates that the CPMG se- 
quence will play a crucial role in extending electron spin 
coherence times in semiconductors. 

In conclusion, we have provided a scheme for coher- 
ence control of localized spins in semiconductors subject 
to general nuclear spin-/ fluctuations. We find that the 
rate for 7r-pulsing needed for substantial coherence en- 
hancement is ^ I^ , yielding an efficient route for achiev- 
ing long time spin coherence in semiconductors using ad- 
vanced pulse technology. 
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